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The circumbinary disk of HD 98800 B: Evidence for disk warping 
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ABSTRACT 

The quadruple young stellar system HD 98800 consists of two spectroscopic 
binary pairs with a circumbinary disk around the B component. Recent work 
by Boden and collaborators using infrared interferometry and radial velocity 
data resulted in a determination of the physical orbit for HD 98800 B. We use 
the resulting inclination of the binary and the measured extinction toward the 
B component stars to constrain the distribution of circumbinary material. Al- 
though a standard optically and geometrically thick disk model can reproduce 
the spectral energy distribution, it can not account for the observed extinction if 
the binary and the disk are co-planar. We next constructed a dynamical model 
to investigate the influence of the A component, which is not in the Ba-Bb orbital 
plane, on the B disk. We find that these interactions have a substantial impact 
on the inclination of the B circumbinary disk with respect to the Ba-Bb orbital 
plane. The resulting warp would be sufficient to place material into the line of 
sight and the non-coplanar disk orientation may also cause the upper layers of 
the disk to intersect the line of sight if the disk is geometrically thick. These 
simulations also support that the dynamics of the Ba-Bb orbit clear the inner 
region to a radius of ~3 AU. We then discuss whether the somewhat unusual 
properties of the HD 98800 B disk are consistent with material remnant from 
the star formation process or with more recent creation by collisions from larger 
bodies. 
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Introduction 



As the majority of stars are formed in binary or higher order multiple systems, under- 
standing the effects of multiplicity is essential to understanding star formation. HD 98800 
is a member o f the TW Hydra assoc iation, whose members have distances from 40 to 100 
pc (Hipparcos; iPerrvman et all Il997l) . and is a quadruple system comprising two spectro- 
scopic binaries ([Torres et al.l 119951 ). HD 98800 A and B, both with K star primaries. The 
AB pair has a s eparation of / . / 8, corresponding to 34 AU and a period of 300-430 years 
( Tokovininl 11999 ) . Ag e estimates for the TW Hya association range from 5 to 20 Myr 
( Soderblom et allll998 , and references therein) and several members including HD 98800 
(IZuckerman fc Becklinlll993l ). have significant infrared excess emission. 



Gehrz et all (1199a ). iKoerner et all (l2000h and IPrato et all (bOOlt ) resolved the AB sys- 



tem excess and determined that the mid-infrared excess arises from material around the B 
spectroscopic binary, the northern component. The excess i s detected only at wavelengths 



longer than 8 microns and is well fit by a 150 K blackbodv (lLow et al 



2000l ; |Prato et alll200ll ). Recent Spitzer IRS observations by lFurlan et all (120071 ) show a small 
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Koerner et al. 



excess starting at 5.5 /im with the majority of the excess at wavelengths of 8 /zm and longer. 
At this temperature the material is several AU f rom the stars and must therefore be cir- 
cumbinary given the semi-major axis of 0.98 AU dBoden et alJbood ). HD 98800 shows no 



sign of active accretion a s judged by the Ha line width of 0.19 Angstroms (IDunkin et al. 
19971 ) or X-ray emission (IKastner et al.l 120041 ) . This is in contrast with two of the other 



association members, TW Hya and Hen 3-600 A, which have observ ed accretion rates onl y 
1-2 orders of magnitude lower than those in 1 Myr-old T Tauri stars (jMuzerolle et al.ll2000l ). 



Circumbinary disks are no t uncommon around clo sely spaced (separation less than a 
few AU) young stellar binaries (I Jensen fc Mathieulll997l ). However, at an age of 10 Myr, 80- 
90% of young stars have d ispersed their primordial disk as traced by near and mid-infrared 
excess ( Hillenbrand 1120071 ). The population of main sequence stars w ith an infrared excess 
(the so-called debris or secondary disk systems) discovered by IRAS (lAumann et al.l 1 19841 ) 
has been greatly expanded by the sensitivity of Spitzer. Recent Spitzer o bservations show 
^ 15% of A to K spectral type field dwarfs have a mid-infrared excess (e.g. iBeichman et al. 



(120061 ); IChen et al.l (120061 )). but the material in these disks is believed to arise from collisions 
between larger bodies rather than being primordial and many of these stars have an excess 
only at wavelengths > 20/im, indicative of cool dust. 

Determining the evolutionary state of the HD 98800 B circumbinary disk is compli- 
cated by multiplicit y. Secondary disks are also observed around main sequence binaries 
( Trilling et al.l 120071 ) as they are around pre-main sequence binaries (jJensen fc Mathieu 
1997h . Dynamical studies have shown that a circumbinary disk will have a cleared cen- 
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tral region which can greatly slow or stop mate rial from accreting onto the central star 
( Artymowicz Sz Lubow 1994 ; Pichardo et all 2005 ). Thus it is not clear if the material sur- 
rounding B is remnant from the star formation process or has been created by collisions from 
larger bodies. 



Recently, iBoden et al.l (120051 ) used measurements from the Keck Interferometer, the 
Hubble Space Telescope Fine Guidance Sensors and radial velocity data to d etermine the 



physic al orbit of HD 98800 B and the mass and luminosity of its components. IBoden et al. 



( 120051 ) modeled the SEDs of the HD 98800 B components, finding effective temperatures 
of 4200 ± 150 K and 4000 ± 150 K, and luminosities of 0.330 ± 0.017 L and 0.167 ± 
0.038 L for the primary and secondary respectively. Additionally, a small ex tinction of 
Ay 0.3 ± 0.05 was required to fit the SED of the system. ISoderblom et all (119981 ) also found 
that a small extinc tion (Ay 0.44 mag ) was necessary to explain the lithium-predicted 
visible photometry. ISoderblom et all (119981 ) additionally noted the significa nt time var i ability 
with no obvious periodicity in the Hipparcos optical photometry, which iTokovininl (119991 ) 
suggests may be due to extinction through an edge-on disk. In contrast, an earlier monitoring 
campaign by iHenry fc Halll (1l994l ) detected a 0.07 magnitude variation at V with a 14 day 
period which they suggest is due to a starspot; however, these authors were unaware of the 
dual spectroscopic binary composition of the HD 98800 system. 

In this paper, we use the orbital parameters and the observed SED to constrain the 
physical properties of the circumstellar material around HD 98800 B. In §31 we investigate 
the origin of the extinction towards HD 98800 B given the inclination of the binary orbit 
using geometrically thick circumstellar disk models. In §H we present a dynamical simulation 
of the interaction between the A and B components and the B disk. In $5j we compare our 
results to other models of the HD 98800 B disk and discuss the implications of the modeling 
on the question of whether the disk material is primordial or secondary and on the origin of 
the extinction along the line of sight to HD 98800 B. Our conclusions are summarized in §3 



2. Input data 



2.1. Photometry 



We use near-infrared t hrough sub-millimeter photometry from the literature (jPrato et al 



200ll ; ISylvester et al.l 119961 ) to determine the flux contribution of the disk. In this work, we 
concentrate only on the continuum emi s sion a nd d o not use the photo metry affected by the 



10 jum silicate feature. iKoerner et al. ( 200 



feature and recent work by iFurlan et all (1200 



and iPrato et all (120011 ) discuss this spectral 
using the Spitzer mid-infrared spectrum sug- 
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gests the presence of amorphous carbon, pyroxene and olivine grains. The stellar effective 
temperatures, luminosities and system distance are taken from iBoden et al.l (120051 ). 



2.2. Orbital parameters 



We use the relevant orbital parameters from lBoden et al.l (120051 ) (Tabled]) to derive the 
dynamical constraints on the disk. Circumstellar material will tend to be cleared from the 
inner disk region through interactions with the central stars (see e.g. lArtymowicz fc Lubow 
19941 ; iPichardo et all 120051 1 and this clearing increases with the orbital eccentricity of the 
binary. For HD 98800 B the apastron distance is 1.7 5 AU, making any m aterial within this 
radius dynamically unstable. From the calculations of lPichardo et all (120051 ) the circumbinary 
disk for a system such as HD 98800 B with eccentricity of 0.78 and secondary mass fraction 
^2/ (wii + — 0.45 is cleared to a radius of 3.5 times the semi-major axis or 3.4 AU. 



e 


0.7849 ± 0.0053 




i (deg) 


66.8 ± 3.2 




a (mas) 


23.3 ± 2.5 




System distance (pc) 


42.2 ± 4.7 






Ba 


Bb 


Mass (Mq) 


0.699 ± 0.064 


0.582 ± 0.051 


T e // (K) 


4200 ± 150 


4000 ± 150 


Luminosity (L Q ) 


0.330 ± 0.075 


0.167 ± 0.038 



Table 1: Orbital and physical parameters for HD 98800 B from IBoden et al.l (120051 ) 



The outer radius was derived using similar constraints from the A-B orbit parameters 
from Tokovinin (1999; a = 62 AU and e = 0.5; Orbit II), an A comp onent total mass of 1.3 

an d the circumsecondary truncation resul ts shown in Figure 7 of lArtymowicz fc Lubow 
(119941 ) for a viscosity of 10~ 3 . Although the lArtymowicz fc Lubowl (119941 ) models are for 
a circumseconda ry disk in a planar orbit and the A-B and Ba-Bb orbits are not coplanar 
(iTokovininl Il999l ) we use this as a guide to how any material distributed around B would 
be dynamically influenced by the A-B orbit and these inner and outer truncation radii are 
confirmed for the case of HD 98800 B by the dynamical model presented in £|4l 
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2.3. Extinction 



Mid-infrared imaging from iKoerner et al.l (120001 ) and iPrato et all (120011 ) indicated the 
large IR excess in HD 98800 is associated with HD 98800 B, with no significant excess 
emission from HD 98800 A. Spectroscopic and interfer ometric obser v ations that resolve the 
flux ratios from the B subsystem components allowed iBoden et al.l (120051 ) to model the B 
subsystem SEP, an d assess the individual component temperatures and luminosities. The 
Boden et al.l (120051 ) SED modeling required a small a mount of extinct i on to match visible 
HST photometry; a similar conclusion was drawn by ISoderblom et all (119981 ) using visible 
spectroscopy. Finally, both groups concluded that no extinction is apparent in the SED of 
HD 98800 A. The association of IR excess and extinction with B and not A suggests that 
there is circumbinary material around the B subsystem. 

To quantify the c olumn density of d ust which would produce Ay 0.3 ± 0.05, we used 
the dust parameters of IWood et al.l (120021 ) derived for the HH 30 disk. These parameters are 
most likely to be appropriate if the extinction arises due to the line of sight passing through 
some portion of the HD 98800 B circumbinary disk and assuming a similar co mposition for 



the H H 30 and HD 98800 B disks. At V (0.55 /mi), the dust mass opacity from lWood et al. 



( 120021 ) is 35 cm 2 /gm. For an extinction of Ay 0.3 ± 0.05, this results in a column density of 
8 x 10~ 3 g/cm 2 . For comparison, the surfa ce density of the minimum mass solar nebula at 
1 AU is ~ 10 3 g/cm 2 (I Weidenschillingl 1 1 9 77l ) . thus if the extinction is due to the line of sight 
passing through the disk, it is only through the tenuous outer layers. 



3. Circumbinary disk models 



Several previo us groups have successfully reproduced the SEP of HP 98800 B with a 
~150 K blackbody JLow et all 19991 : IKoerner et all200olJPrato et alllioOll : iFurlan eTa1l2007h . 
However, none of these groups includ ed all the dynamic al constraints (3 < R s tabie ^10 AU) 
and the unresolved mid-IR imaging of IPrato et all (120011 ). Here we consider if a geometrically 
thick dust and gas disk can reproduce the SEP while fitting both the dynamical constraints 
and intersecting the line of sight to provide the observed extinction. 

We use a Bayesian approach to characterize the disk pa r amete rs; t his approach has bee n 
used to characterize disks around young stars by lLay et al.l (119971 ) and lAkeson et al.l (120021 ). 
We fix some of the physical characteristics based on the assumption of the circumbinary disk 
being co-planar with the B binary orbit (i.e. inclination) and other parameters based on the 
physics of the input model (i.e. radial temperature gradient). For the remaining parameters 
we use the Bayesian formalism to calculate a probability distribution for each parameter. 
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This method shows how well a given parameter is constrained and how the free parameters 
may be degenerate with each other, which is a particular concern using only spatially unre- 
solved data. It also facilitates the inclusion of additional observational constraints (e.g. the 
dynamical truncation radii). In our case, the assumptions made (the inclination angle and 
co-planarity of the disk and the Ba-Bb orbit) do not have associated probabilities, so the 
relative likelihood a single model, given the known data set, is proportional to e - * 2 / 2 . This 
probability is calculated for every possible set of parameter values. Then, for each parameter 
value, the relative probability is the sum of all probabilities for the range of values in all the 
other parameters. 



3.1. Flared accretion disk 



We first consider a disk with physical properties relevant to the active accretion phase 
of younger T Tauri stars and in particular, consider a flared disk model to investigate if the 
extinction can arise from the line-of-sight intersecting the outer regions of a flared disk. A 
simple calculation shows that this is plausible. The exact form of the disk height (H) as a 
function o f radius (r) depends on the disk temperature profile a nd ranges from H oc r 9 / 8 to 
H oc r 9//7 ( Kenvon fc Hartmannlll987 : IChiang fc Goldreichlll997l ). Here, we use the relations 
derived by lchiang &: Goldreich ( 1997 ) where at radii relevant to the HD 98800 B disk H = 
0.l7r 9 ' 7 where r is the disk radius in AU. Theref ore the angle of inci dence (a) for light from 
the central star is tana = H/r = 0.17r 2 / 7 . From Boden et al. ( 20051 ). the orbital inclination 
angle is 67°± 3° and therefore a > 23° when r > 25 AU (Figured]). Although an outer radius 
this large is probably ruled out by the A-B orbit truncation effects, the observed column is 
so small that it may arise from material above the disk scale height. 

To further investigate the possibi lity that the extinction ari ses from the outer regions 
of a flared disk we parameterized the IChiang fc Goldreich! (119971 ) disk and determined the 
disk properties as constrained by the SED data. The disk is comprised of an interior and a 
superheated surface. For the wavelengths and radii under consideration in the HD 98800 B 
disk, the surface layer is always optically thin and so the surface flux is simply added to the 
interior flux (which may be optically thick). To account for a different stellar flux than in 
the original model, we use the interior Tj and surface T s temperature radial scalings from 
Chiang fc Goldreich! (119971 ). but include an single scaling factor T , 

Ti = T o (150/r 3 / 7 ) K (1) 

T s = T o (550/r 2 / 5 ) K, (2) 

where r is the radial distance in AU. The surface density (£) and its radial exponent (p) are 
also allowed to vary. The resulting model has six parameters: T Q , S G (surface density at 2 
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Fig. 1. — Illustration of the angles and scales for the HD 98800B system. The figure is 
looking edge-on to the Ba-Bb orbit (thick solid line) with the center of mass at the origin. 
The inferred line-of-sight to Earth is shown as a dotted line and the ang le with the A-B orbit 
i s show n as a dashed line. The scale height for a flared disk following the lChiang fc Goldreich 
(119971 ) calculations ( §3.11) is also shown, with an inner radius of 3.4 AU. The outer radius is 
deliberately unconstrained to demonstrate where the disk would intersect the line-of-sight if 
the radial extent were sufficiently large. 



AU), p, r in , Ar (disk radial extent) and (3 (dust opacity frequency exponent) and the range 
of values searched for each is given in Table [2j The inclination angle is set at 67°. 

All parameters have a range of values which are contained within models which roughly 
match the SED. Table 2 gives the 90% probability range for each parameter. We note that 
these models tend to produce too much flux in the mid-infrared, which can be eliminated 
by removing the surface layer and having a single layer disk; however, the probability dis- 
tribution for the inner radius for the single layer models does not include any models where 
r in > 3 AU (fitting the dynamical constraints). It may be possible to produce a better fit to 
the mid-infrared SED and satisfy the dynamical constraints with modifications to the tem- 
perature structure of the disk ; but we did not wish to add mor e parameters. In comparison 
to the single temperature fits (ILow et al.lll999l ; iPrato et alll200ll ) the interior temperatures at 
the inner radius are below 150 K, but this is balanced by the much hotter surf ace layer. An 



example model, along with the SED data used for the stellar photosphere fit (IBoden et al. 
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Parameter Range searched 90% prob. range 



T 


0.35 - 


-0.7 


0.45 - 0.6 


£ (gm/cm 2 ) 


0.02 


- 3 


0.025 - 0.5 


P 


0.5 - 


2.5 


unconstrained 


Ar (AU) 


0.2 


- 4 


0.8 - 2.0 


r m (AU) 


1 - 


4 


1.1 - 3.6 


P 


- 


2 


0.15 - 0.9 



Table 2: Range of values searched for each parameter of the flared disk model and the range 
of values which encompasses a probability range of 90%. 



(120051 )) and the data used in the disk modeling, are shown in Figure [2J The SED can be 
roughly reproduced with a range of inner radii from 1.1 to 3.6 AU. Interestingly, all the 
probable models have a relatively small radial extent (Ar < 2 AU) with the 90% proba- 
bility range from 0.8 to 2.0 AU. A d isk with this extent would appear un resolved in the 
mid-infrared imaging (consistent with lKoerner et al.ll2000t iPrato et all 120011 ) . 



We can now calculate if the models which reproduce the SED can also explain the 
observed extinction. Using an example model with rj n = 3.2 AU, Ar = 0.8 AU and a 
surface density of 0.035 gm/cm 2 at a radial distance of 2 AU, the column density along 
the derived line of sight is ~ 10~ 12 gm/cm 2 , with similar column densi t ies de rived for other 
models with high relative probabilities. For the IChiang fc Goldreichl ( 119971 ) model visible 



dust opacity of Ky ~ 400 cm /gm, this results in an extinction of Ay ~ 10 , many orders 
of magnitude below the observed extinction of 0.3 ± 0.05. 



3.2. Vertically extended inner radius 



Next, we considered if the extinction could arise from a vertically extended rim at 
the inner disk rad i us. Vertically-extended inner rims arise from direct ir radiation of dust 



( Natta et al.ll200ll ; iDullemond et al.l l200ll ) at the inner disk edge. From iDullemond et al. 
( I200ll ) the pressure scale height is 



JcT ■ R 3 

rv - L rim 1 L rim 
VgGM* 



1/2 



(3) 



where h rim ,T rim and R r 



are the pressure scale height, temperature and radius of the inner 
rim, fig is the mean molecular weight of the gas and the rim height H r i m 5h r i m . Using a 
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10 100 
Wavelength (/xm) 



1000 



Fig. 2. — The stellar photosphere models for the components of HD 98800 B (dashed = 
Ba, dotted = Bb, solid=total) along with photometry from the literature. The optical to 
near-infrared photometry was used for the stellar fit, while the mid-infr ared to mill i meter 
pho tometry was use d for th e disk modeling. The photometry is taken from lPrato et al ( 2001 ) 
and ISylvester et al.l (119961 ) and where error bars are not visible, they are smaller than the 
symbol size. An example flared disk model (dot-dashed) is shown with T o =0.61, S = 0.035 
gm/cm 2 , p = 1, r in = 3.1 AU, Ar = 0.8 AU and (3 = 0.8. 



combined mass of 1.28 M , \i g = 2 and T rim =150 K, 



±J -rim ~ 



0.12 



R, 



AU 



1/2 



AU. 



(4) 



If the inner disk is dynamically truncated by the binary orbit at a radius of 3.4 AU ( § 2 . 2 [) . 
the height of the rim is 0.41 AU. A rim at this radius and height would obscure the line of 
sight starting at an inclination angle of 85°. Given the inclination angle of the binary orbit 
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of 67° ± 3°, the inner disk rim will not obscure the line-of-sight if the disk is co-planar with 
the binary orbit. 



3.3. Circumbinary envelope 

Finally, we consider the possibility that the extinction towards HD 98800 B could arise 
from material distributed in a circumbinary envelope around the B pair. Assuming a spheri- 
cal distribution, the surface density calculated above can be used to estimate the total mass 
in such an envelope if the inner and outer radius are known. For the purpose of a rough 
mass estimate, we use a value for the inner radius from the region cleared by the B orbit of 
3.4 AU and a value for the outer radius of 10 AU ( §2.21) . which results in an envelope mass 
©■ 



of 5 x KT 7 M 



Circumstellar envelopes are common around young stars wi th ages of a few Myr but 



have generally dissipated at the derived age of HD 98800 B (e.g., iMundy. Looney. fc Welch 



20001 ). However, the dynamics of the quadruple system are significantly more complicated 
than considered in the canonical star formation picture and a small amount of material from 
the formation of the central stars may have been unable to accrete onto the disk. We consider 
this unlikely given that the A component, which has a similar mass and is presumably coeval 
has no such material. 



4. Dynamical simulations 

The HD 98800 B circumbinary disk exists within the dynamical environment of the 
larger quadruple system and could be warped by interactions with HD 98800 A. To model 
the HD 98800 dynamical environment, we consider a system comprising 3 stars. Two of these 
stars represent the HD 98800 B system with masses of M x = O.699M and M 2 = O.582M , 
an orbital period of 315 days and an eccentricity of e = 0.78. These two stars are then 
surrounded by a circumbinary disk comprising 25000 test particles with an initial inner 
radius of 2 AU and an outer radius of 15 AU. The disk particles were distributed randomly 
with no initial eccentricity, but with initial inclinations that were Rayleigh-distributed, with 
the sine of the inclinations, sinz, having a root mean square value of 0.05. The disk is 
therefore initially thin, but not completely flat. 

The HD 98800 A binary system was then represented by a single star with a mass equal 
to the combined mass of the two stars representing the HD 98800 B system (M 3 = 1.281M ). 
Although we could have modeled HD 98800 A using two stars, the actual orbital parameters 
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and stellar masses are not well known, and this would have been computationally expensive. 
Representing HD 98800 A with a single star is an approximation that should at least give 
us some idea of how this outer system influences the disk around HD 98800 B. The orbital 
perio d of the A-B system was taken to be 345 years with an eccentricity of e = 0.5 (ITokovinin 
19991 . Orbit II) and with an orbital plane initially inclined at 10° to the plane of the disk 
around HD 98800 B, which is coincident with HD 98800 B's orbital plane. 

I n this dynamical si mulation we have ignored radiation pressure and Poynting- Robertson 
drag (IBurns et al.lll979l ). HD 98800 is a young sy stem and the disk around HD 98800 B is 
relatively massive and is collisionally dominated (jWyatt et al.ll2006l ). Poynting- Robertson 
drag in a system like HD 98800 B is therefore expected to be insignificant because dust 
will be des troyed by col lisions before Poynting-Robertson drag can cause significant inward 
migration jWyatt l l2005h . Altho ugh radiation press ure could change the orbital parameters of 
particles produced by collisions (jWyatt et al.lll999l ). the pericenters of these orbits should be 
randomly orientated and so the effect will be to spread these particles out without changing 
the plane of their orbits. Since this is unlikely to affect significantly the global structure of 
the disk we ignore it here and include only the gravitational influence of the stars. 



The system was evolved for just over 10 6 years using a Hermite integrator (IMakino 



199ll ). that conserved both energy and angular momentum to within 0.1%. We found that 
the star representing HD 98800 A caused the HD 98800 B binary system and its circumstellar 
disk to undergo torque induced precession with a period of ~ 99000 years. The inclination 
of the HD 98800 B system, with respect to the initial plane of the HD 98800 A/B orbit, also 
varies sinusoidally from the initial 10° to just over 20° with a period of ~43000 years. 

Figure [3] shows the variation of the disk structure over one precession period from 
t = 1.0823 Myr to t = 1.1841 Myr. In each panel the center of mass of the HD 98800 B 
system is located at the origin, and each panel shows a thin slice through the circumbinary 
disk together with the central stars of the HD 98800 B system. In these figures, the x-y 
plane is the initial plane of the orbit of the HD 98800 A/B system. The orbital plane of 
HD 98800 B and its circumbinary disk were therefore initially inclined at 10° to the x-y 
plane. The label on each figure shows the time, the viewing angle ((f)), and the inclination of 
the orbital plane of HD 98800 B relative to the x-y plane (9). The viewing angle is varied to 
take into account the precession of HD 98800 B, such that the disk is always viewed along 
the axis, in the x-y plane, about which the disk appears to have been rotated. 

What the panels in Figure [3] show is that the eccentric orbit of HD 98800 B (e = 0.781) 
truncates the inner disk at a radius of ~ 3 AU (in agreement with the estimates in §2.2p . 
while the outer disk, which initially extended to a radius of 15 AU, is now truncated at a 
radius of ~ 10 AU due to the influence of HD 98800 A. Note that this outer truncation radius 
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-10 -5 5 10 -10 -5 5 10 -10 -5 5 10 

si (All) x (AU} x (AU) 

Fig. 3. — A series of figures showing the evolution of the circumbinary disk around 
HD 98800 B over one precession period. The center of mass of HD 98800 B is located 
at the origin and the x-y plane is the initial plane of the HD 98800 A/B orbit. Each panel 
shows a slice through the circumbinary disk together with the two central stars. The viewing 
angle (0) is varied to take into account the precession of the system and such that the disk 
is always viewed along the axis, the in x-y plane, about which the disk appears to have been 
rotated. Not only do these figures illustrate how the inner and outer disk are truncated, they 
also illustrate how the disk becomes warped as the angle, 9, between the orbital plane and 
the x-y axis increases. This warping is qualitatively consistent with the visual extinction 
observed towards HD 98800 B, which has an inclination angle of 67°. 

also agrees with the estimates based on co-planar interactions ( §2.21) . This outer radius value 
further supports the idea that the observed visual extinction to HD 98800 B is unlikely to 
be due to a flared co-planar disk extending to large radii. What is also clear from Figure 
[3] is that when the inclination of the orbital plane relative to x-y plane is small (6 ~ 10°), 
the disk is thin and slightly flared. However, as the angle increases, the disk becomes more 
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obviously warped. Consider the third panel in the top row. The orbital plane is inclined at 
21° to the x-y axis and yet a line of sight along the x-axis would still intercept some disk 
material. 

Although the resolution of our simulation is insufficient to quantify the column density 
of material located more than 20° above the orbital plane of HD 98800 B, the existence of 
the warp is qualitatively consistent with the observed extinction to HD 98800 B, the orbit 
of which has an inclination angle of 67°, particularly given the low column density required 
to produce the observed extinction. T his is consistent with the observed variability in the 
Hipparcos data (jSoderblom et allll998l ) arising from changes in extinction if the disk material 
is clumpy, but as these simulations can not quantify the column density, we can not directly 
compare to the variability. The rot ation time scale of th e warp is much longer than the 14 
day periodic variations observed by iHenry &: Hall! (119941 ) which does not correspond well to 
the disk rotation time scale of 2.5 yrs at 3 AU. 



5. Discussion 

In this section, we first compare our results on the circumbinary disk of HD 98800 B 
to models of this disk from other authors. Next we compare the HD 98800 B disk to 
samples of other objects with primordial and secondary disks, discuss the implications for 
the evolutionary status of this object and suggest some observational tests of this status. 



5.1. Comparison to other models 



As listed in jjj and jj3l several other groups have modeled the dis k of HD 98800 B 
JLow et alill999l : iKoerner et aDl200ol : IPrato et allboOll ; iFurlan et ajl2007h . All these models 
have the following characteristics in common: a circumbinary disk with a inner region cleared 
to at least a few AU do minated by emissi on fr om material with a temperature of ~150 K. 
Here we compare to the IPrato et all (120011 ) and IFurlan et all (120071 ) models which use either 
imaging constr aints or deta i led sp ectroscopy in their modeling. Using their unresolved images 



as constraints, IPrato et all (120011 ) fit the SED with a blackbody temperature of 150 K and 
derived an inner radius 2 AU, an outer radius of at least 5 AU, and a possible cooler dust 
component at larger radii which is then truncated by the A/B component interactions. They 
also comment that the vertical extent must be ~1.5 AU to account for the hig h fract ional 
luminosity. Our dyn amical simulations and the work of lArtymowicz &: Lubowl (119941 ) and 
Pichardo et all (120051 ) suggest that this inner radius is too small to be stable given the orbital 
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parameters of the binary. 



Recently, iFurlan et all (120071 ) used the mid- infrared spectrum of HD 98800 from the 
Spitzer telescope to constrain the structure and composition of the HD 98800 B disk. Their 
spectrum shows no excess below 5.5 fim (within a 5% absolute calibration) and a strong 
excess starting at 8 /xm. They use the structure in the silicate feature to model the dust 
composition and we refer interested readers to the paper for more details on that aspect. 
In comparison to our model, the relevant model parameters are an optically thick region 
starting at 5.9 AU with a temperature of 150 K and an optically thin ring between 1.5 and 2 
AU with grain sizes of a few microns and a temperature of 310 K where the silicate emission 
features arise from. This inner ring is within the apastron separatio n of the Ba-Bb pair of 
1.75 AU and within the dyn amically cleared regions in the work of lArtymowicz fc Lubow 
(I1994J ). iPichardo et all (120051 ) and §H If material is located at these radii, it will be transient 
as the time scales for clearing the gap (<100 orbital periods, i.e. <100 yrs) are much shorter 
than the P-R drag timescale of 10 5 yrs for a 3 /zm grain at 3 AU in this system. HD 98800 B 
also shows no indication of accretion, unlike the DQ Tau and UZ Tau E syst ems which are 



inferred to have active ac cretion through streamers from a circumbinary disk (iMathieu et al. 
19971 : 1 Jensen et al.ll2007h . 



Both lPrato et all (120011 ) and lFurlan et all (120071 ) noted that small grains must be present 
in the HD 98800 B disk to account for the spectral features. A distr ibution of grains s ize is 
inferred in many pre- main sequence, primordial dis ks (see review bvlNatta et al.ll2007l ) and 
the higher temperatures silicate grains inferred by IFurlan et all (120071 ) from features in the 
mid-infrared spectrum could reside on the stellar photon-illuminated surface of a flared or 
warped disk such as that discussed in §3.11 In the example model presented in Figure 2, the 
interior temperature at the inner rim is 60 K, while the surface layer temperature is 210 K. 
Spectrally resolved mid-infrared interferometry observations would constrain the location of 
the small grains responsible for the silicate emission feature. 



5.2. Comparison to other primordial and secondary disks 



The TW Hya association is unusual in containing sources with a very large range of 
pre-main sequence properties. Four of the sources have JiiMS-detected infrared excesses 
(TW Hya, HD 98800, Hen 3-600 and HR 4796A). TW Hya and Hen 3-60 are often classi- 
fied as T Tauri stars given their accretion signatures dMuzerolle et al.ll2000l) and near-infrared 



exces s, and TW Hya has observed CO and H2 in the disk (IDent et al.ll2005l ; IWeintraub et al. 



2000( 1 . On the other hand, HR 4796A has a similar SED to HD 98800, but with material 
at much larger radii and is consistent with a secondary dust debris disk. Thus, it is not 
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clear from the association as a whole whether the HD 98800 B disk is primordial or sec- 
ondary. Here we compare the derived disk properties of HD 98800 B with those of other 
primordial/protostellar and secondary/debris disk systems. 

The mass distribution of the most probable models from §3.11 for a flared gas and dust 
disk peaks at 10 -7 M , significantly lower than most classical T Tauri stars an d comparable 



to ma ny secondary debris disk stars such as the beta Pic moving group (e.g. iMeyer et al. 
(120071 )); however, as shown in §U the HD 98800 B circumbinary disk is likely to be dynam 



ically truncated at both the inner and outer radii, and the disk mass is highly correlated 
with the small radial extent (few AU) as compared to values of ~100 AU typical for T Tauri 
circumstellar disks. An additional contribution to the low disk mass is the surface density 
of <0.5 gm/cm 2 at a radius of 2 AU, which is a factor of ~100 lower than the minimum 
mass solar nebul a. The high infrared luminosity of the disk (0.1 - 0.2 L^) has been noted 



by many groups (IZuckerman fc Becklinlll993l ; iKoerner et al.ll2000i ; iPrato et alll200ll ) and this 
fra ction is significantl y higher other young debris disks systems, such as the sample observed 
by IChen et al.l (120061 ) of A to M spectral types with ages <50 Myr. In that sample, the 
detected systems had infrared disk luminosities of 5 x 10~ 6 — 5 x 10 -3 L/^/L*. However, the 
dust in those systems is in general located many to tens of AU from the central star, rather 
than a few AU and is not likely to be as vertically extended. Overall, the disk properties for 
HD 98800 B are intermediate between those of T Tauri, primordial disks (and its association 
members TW Hya and Hen 3-600) and secondary debris disks around stars with ages of tens 
of Myr. 

A more stringent diagnostic of primordial vs. secondary material would be to observe 



gas tr acers in the disk and establish the gas to dust ratio. In a CO survey by iDent et al. 



(120051 ). HD 98800 was not detected with an upper limit of 0.09 K km/sec, while TW Hya was 
detected with 2.16 K km/sec of emission. Although their relative infrared excesses (~0.1 
for HD 98800 and 0.3 for TW Hya) are similar, the difference in CO emission strengths 
may result primarily from the CO arising in colder gas 10's of AU from the star. As the 
HD 98800 B disk is probably much smaller than this, the CO non-detection does not set 
stringent limits on gas depletion. 

The lack of accretion signature in HD 98800 B is not surprising in the context of a 
~3 AU inner hole; however, accretion across the gap in a circumbinary disk can take place 
(lArtymowicz fc Lubowl 119961 ) , particularly for geometrically thick dis ks and has been ob- 
served in T Tauri circumbinary disks such as DQ Tau and UZ Tau E dMathieu et al.lll997l ; 
Jensen et al.l 120071 ). Dynamical modeling by lArtymowicz fc Lubowl ( 119961 ) shows that the 
mass accretion from circumbinary disks is strongly modulated in time and peaks at perias- 
tron. Thus another observational test for the HD 98800 B disk would be time monitoring of 
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accretion diagnostics, particularly centered about periastron. 



6. Summary 

HD 98800 B is one of the few young multiple stars with significant disk emission for which 
the physical orbit has been determined. We have used the orbital parameters, particularly 
the inclination, and the assumption of disk-bina r y cop lanarity to examine physical models 



for the disk. As first discussed by iPrato et all (120011 ) the large mid-infrared excess and 



unresolved images of HD 98800 B seem to require an optically and geometrically thick disk. 
Either a flared disk or a vertically extended inner rim can satisfy the geometric requirement, 
but dynamical truncation will keep the majority of the material beyond an inner rim with 
a radius of 3 AU. The disk solution presented here is not unique in reproducing the SED 
and somewhat over-produces mid-infrared continuum flux, but does satisfy the dynamical 
constraints. 

We find that a standard disk prescription is unable to explain the observed extinction 
if the disk and the binary orbit are co-planar. For the flared disk model in Figure [2] to 
intercept enough of the line of sight to produce the observed extinction, the disk inclination 
would need to be 81° as compared to the binary orbit inclination of 67°. From the dynamical 
simulations shown in §H interactions between the B circumbinary disk and the A component 
are sufficient to perturb the disk inclination by that amount. They are also sufficient to warp 
a thin disk enough to place material into the observed line of sight. In either case, we suggest 
that the small observed extinction is due to significant perturbations of the circumbinary disk 
by the A binary component in this quadruple system. The next step in understanding the 
HD 98800 B disk would be to quantitatively combine the SED and observational imaging 
constraints (unresolved mid-infrared size and extinction along the line of sight) with the 
dynamical interactions between all the components in the HD 98800 system. 



We thank Sverre Aarseth for help with the use of his N-body codes and Kenneth Wood 
for initial SED modeling. We also thank Peter Brand and Mark Wyatt for useful discussions 
about secular perturbation theory. We are grateful to the anonymous referee for several 
helpful suggestions. 
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